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ABSTRACT 


Cirrus forecasting is of particular importance to various Department of Defense 
programs. This thesis takes a case study approach to study Air Force Weather Agency 
Mesoscale Model 5 (AFWA MM5) skill in forecasting cirrus clouds, which are not 
represented explicitly by the model (ice water mixing ratio is used as a surrogate.) Two 
cases are selected for study. For each case, an initial forecast time of interest is 
determined which serves as the beginning point for the case study. GOES data and 3- 
hourly MM5 data are then obtained at 3-hourly intervals to coincide with model forecast 
time steps between the initial time through the 30-hour forecast. A standard analysis is 
performed on all data to determine general atmospheric structure for each case at each 3- 
hourly point. Following this, the model’s relative humidity with respect to ice, explicit 
ice water content, vertical velocity, and other fields are considered to determine if the 
model possesses the proper dynamical factors for cirrus formation. Finally, model 
coverage of ice cloud is compared to the ABE cloud mask results to determine how well 
the model’s ice cloud forecasts verify against each 3-hourly observed ice water field 
taken from the GOES data. 

Results indicate that the MM5 underforecasts cirrus coverage, and that the 90% 
relative humidity field with respect to ice may be a better approximation of observed 
cirrus coverage than the ice water field. 
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I. INTRODUCTION 


Mesoscale models have long been in use for operational weather forecasting, yet 
verification of cirrus cloud forecasts by such models has received relatively little 
attention compared to verification of other sensible weather phenomena near the surface 
such as wind gusts, coastal jets, and precipitation. There are several reasons this is the 
case: verification of weather near the surface is relatively easy due to the higher number 
of surface weather reporting sites; detection of surface-based features is relatively 
inexpensive; and sensible weather at the surface has an effect on a wider number of 
people and missions and is therefore of greater interest than verification of weather 
effects aloft. Additionally, as more observing sites become automated, the inability of 
automated stations to report clouds above 12,000 feet results in fewer stations reporting 
high clouds. When mission-impacting weather is considered, we generally do not thi nk 
of cirrus, because cirrus clouds do not adversely affect most day-to-day activities. 

Cirrus’ effects are not often considered important for most of our routine 
activities, but cirrus can have a profound effect on laser propagation. Understanding 
these impacts, and knowing where cirrus will be and in what amounts, is critical in 
providing accurate support to various Department of Defense (DoD) programs that rely 
on lasers. Such programs include the Transformational Communications Architecture 
(TCA) and the Airborne Laser (ABL) Program, as well as other projects of the High 
Energy Laser Joint Technology Office. 

The uncertain nature of basic weather forecasting further complicates forecasting 
the development, presence, and movement of clouds. Many variables affect the 
development and movement of clouds (available moisture, dynamics, microphysics, etc.) 
and there are almost infinite combinations of these variables that can affect clouds in a 
model. A given amount of moisture in an area may assist development of cirrus in that 
area on one day, but may be insufficient for cirrus formation on a different day with 
different atmospheric dynamics. Importantly, many dynamical and microphysical 
processes involved in cloud formation are not fully represented in models. Such 
processes occur on a scale that is well below the model’s ability to explicitly resolve 
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them; therefore parameterizations are developed that attempt to aoeount for these 
proeesses. Parameterizations, however, eannot aeeount for all the details neeessary to 
provide aecurate foreeasts due to their laek of proper representation of mierophysieal 
proeesses. 

Despite the ehallenges posed to eloud foreeasting, MM5 model results are used in 
notional mission planning foreeasts for the programs mentioned above, and similar model 
foreeasts will be of eontinued importanee for those programs as well as future laser- 
related operations. The Air Foree Researeh Lab proposed this thesis research topic to 
address the issue of mesoscale model verification of cirrus cloud forecasts. 

A, BASIS FOR THIS STUDY 

Although not directly related to the meteorological details of this thesis, Ou 
(2002) provides some of the background as to why accurate model predictions of cirrus 
clouds are important when laser-related operations are considered. Ou (2002) developed 
a model to study the effect of backscattering on laser propagation through cirrus clouds, 
with an ABL-centric mindset. He found that laser transmission uncertainty due to cirrus 
particle size and orientation is negligible, but optical depth of cirrus coverage can greatly 
affect laser transmission uncertainty. From observation of cirrus, we know that cirrus 
clouds are not homogeneous with respect to optical depth. Ou (2002) tested two cases 
using homogeneous clouds before comparing those results to a scenario that involved 
real-world satellite data representative of inhomogeneity within cirrus. 

The first model run was done for an aircraft above cloud level with backscattered 
and target-reflected power as functions of target height. Ou (2002) found that 
backscattered power is virtually constant near the aircraft, regardless of cloud optical 
depth, and that target-reflected power is initially much less than the power returned via 
backscattering. As the target left the surface, reflection of laser energy incident upon the 
target increased due to decreasing absorption with height by water vapor in the 
atmosphere. Target reflection decreased as the target approached and entered the clouds 
due to increasing optical depth. Once the target left the top of the cloud, reflected energy 
from the target increased to its maximum value. 
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For an aircraft within the clouds, backscattering was dominated by the 
contributions to backscattering made by the cloud ice particles near the aircraft. Ou 
(2002) found that backscattered energy was one to two orders of magnitude larger than 
for the aircraft flying outside of clouds. The effects of absorption, due to cloud particle 
scattering, aerosol presence, and air molecules, were small. Ice crystal shape and 
distribution greatly affected backscattered power. When the target was below cloud 
level, reflected energy was miniscule when cloud optical depth was 0.2 and was still quite 
small when below a cloud with optical depth of only 0.05. Approaching, within, and 
departing the cloud, the target reflection reached a minimum. Once above the cloud, the 
target reflected more and more energy as atmospheric water vapor decreased with 
altitude. 

After running his model using homogeneous cases, Ou (2002) ran a comparison 
simulation for the case of inhomogeneous cirrus using observed satellite data to deduce 
mean effective ice crystal sizes from a cirrus case observed in northern Oklahoma on 
18 April 1997. Ou placed the aircraft at 11.2 km altitude, 0.1km above cloud top. The 
cloud base was at 7.5 km, and the target began its ascent from the surface. As with the 
homogeneous cases, reflectance of the target increased as it departed the surface, up to a 
maximum sub-cloud value approximately halfway between the surface and the cloud 
base. Values of direct transmission below 4 km were slightly larger for the 
inhomogeneous cirrus event than for the modeled homogeneous event due to the along- 
beam optical depth being shorter for the inhomogeneous cloud. Between 4 km and cloud 
base level, direct transmission was slightly larger for the homogeneous case because of 
the smaller laser beam path through the homogeneous cloud. Although differences in 
transmittance and backscattering were relatively small outside of clouds, scattering and 
direct transmission were both significantly larger when the aircraft was located within the 
modeled homogeneous cirrus region due to the larger laser beam path as well as the 
smaller ice crystal size compared to the case of inhomogeneous cirrus. Below cloud 
level, inhomogeneous cirrus had a very small effect on laser transmission and 
backscattering. 

Ou’s (2002) results show that, for cirrus observed under realistic conditions, a 

missed cirrus forecast can result in ABL mission failure. The presence of unforecast 
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cirrus can significantly reduce laser transmission between aireraft and target, rendering 
the ABL useless in thwarting a threat, or leaving other laser systems inoperable. 

B, OBJECTIVES 

This study seeks to evaluate the AFWA MM5’s ability to accurately foreeast 
cirrus coverage over a 2-dimensional area from a starting time of interest up to the 30- 
hour foreeast period. The study endeavors to obtain a basic understanding of model 
strengths and weaknesses in forecasting ice clouds for weather events a forecaster would 
experienee on a routine basis. 

Speeific goals of this study are to: 

1. Determine whether cirrus forecast errors are mostly due to basie weather 
system forecast errors (e.g., timing or phase errors) or eloud process errors; 

2. Assess how model dynamies, moisture distribution, and eloud microphysics 
contribute to cirrus formation in the model; 

3. Identify possible limitations in MM5 cloud parameterizations that impact on 
cirrus forecast accuracy; 

4. Identify eritieal parameters that might help prediet eirrus even when explieit 
model foreeasts do not report ice water for a given area. 
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II. BACKGROUND/LITERATURE REVIEW 


A. CHARACTERISTICS OF CIRRUS 

It is important to define what is meant by eirrus, what synoptie features typieally 
drive eirrus formation, and what eonditions (not synoptically-driven) ean result in the 
formation of eirrus. 

1, Composition 

Cirriform elouds are eomposed almost entirely of iee erystals that form from the 
freezing of supereooled water droplets. Three genuses of high clouds exist: cirrus 
(cirriform appearance), cirrostratus (stratiform appearance), and cirrocumulus 
(cumuliform appearance). 

Cirrus clouds are commonly known as “mares’ tails”, and appear as strands or 
filaments of clouds due to advection of ice crystals downstream by high winds aloft. 
Stratiform and cumuliform clouds appear as patches; cirrostratus has a smooth 
appearance, while cirrocumulus has a bumpy appearance. Within each genus, several 
species and varieties exist that further define cloud characteristics. This study will use 
the term, “cirrus” to generally refer to ice clouds, regardless of genus. 

2, Location and Related Synoptic-Scale Features 

In mid-latitudes, cirriform clouds are normally found between 20,000 to 30,000 
feet. It is not uncommon to find cirriform clouds below 10,000 feet in polar areas and as 
high as 60,000 feet or more in the tropics. 

Cirriform clouds are typically associated with various synoptic conditions: in 
advance of warm fronts, “spilling over” ridges, associated with jet streams, and lee of 
mountain ranges. Related dynamical details as applicable to the case studies in this 
document are provided in the next section. 

3, Formation Conditions and Basic Governing Dynamics 

Development of cirriform clouds first requires cooling a parcel of air to 

saturation, usually by gentle lifting along an isentropic surface. Following saturation, the 

nucleation phase occurs by one of two mechanisms: heterogeneous nucleation or 

homogeneous nucleation. Heterogeneous nucleation requires a catalyst or ice nucleus 

and typically occurs when temperatures range between -5 degrees Celsius to 

5 



approximately -40 degrees Celsius depending on the availability of iee nuclei. Ice nuclei 
generally take the form of dust particles or preexisting ice particles. Homogeneous 
nucleation occurs at temperatures near -40 degrees Celsius and does not require the 
preexisting presence of ice nuclei. According to Houze (1993), cirrus can exist down to - 
85 degrees Celsius. Typical vertical velocity values for the formation of cirrus are on the 
order of 0.1 to 0.2 m s"' in most cases, although formation of certain convective species 
of cirrus is associated with vertical velocity values of 1 to 2 m s‘\ During the FIRE 
cirrus experiments, Starr and Lare (1993) observed cirrus in areas of vertical velocity 
values on the order of 0.02 m s'\ 

Cirrus clouds are typically located ahead of surface warm fronts where warmer air 
at low levels rises, cools adiabatically, and eventually condenses. Following saturation, 
continued cooling leads to the formation of ice crystals via one of the two methods 
presented above. 

This thesis will focus on two cases of interest. In the first case the cirrus initially 
resulted from the lifting of moisture by convection in the tropics, which then streamed 
into northern Mexico and the desert southwest by the forcing of an upper-level ridge to 
the east and an upper-level closed low off the United States west coast. In the second 
case, cirrus were formed by lifting ahead of an upper level low. These cirrus clouds 
moved across an upper-level ridge downstream of their initial moisture source. In lieu of 
a method to directly observe model parameterization details as time progresses, this 
thesis will investigate the dynamics responsible for cirrus formation in both cases and 
attempt to draw conclusions regarding how well the MM5 places cirrus related to 
modeled dynamical forcing. 

4, Microphysical Properties 

As described by Houze (1993), cirriform clouds typically possess ice contents on 
the order of 0.001 to 0.25 g m' , with the most typical values on the order of 0.01 to 
0.1 g m' . Particle sizes range from 50 to 1000 pm. A variety of shapes can describe the 
ice crystals that compose cirrus clouds, including columns, bullets, bullet-rosettes, and 
plates. Aggregates of the various crystal shapes are fairly common, since partially melted 
ice crystals do occasionally re-freeze upon making contact with frozen crystals (i.e., 
heterogeneous nucleation.) 
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B, CLOUD MASK THEORY 

The cloud mask used by the ABL was developed by Atmospheric and 
Environmental Research Incorporated (AER). The algorithm makes use of a variety of 
tests to determine cloud coverage and phase using various data sources. These data 
sources include (among others) the Air Force Weather Agency Mesoscale Model 5 
(AFWA MM5), the Geostationary Operational Environmental Satellite (GOES), and the 
Moderate Resolution Imaging Spectrometer (MODIS) satellite. In addition to cloud tests, 
background tests are also conducted to determine the presence of sun glint, desert 
backgrounds, and snow/ice backgrounds. Knowledge of background conditions is critical 
to successful results from certain cloud tests. 

Although data from the MODIS hardware can also be used, the cloud mask for 
this thesis was run strictly on GOES satellite data (specifically, GOES-10.) The GOES- 
10 imager contains sensors centered at the 0.65, 3.9, 6.7, 10.8, and 12.0 pm wavelengths. 
The GOES-12 imager also has the 0.65, 3.9, 6.7, and 10.8 pm wavelengths in the first 
four channels, however the 12.0 pm band is no longer present, and a 13.0 pm band has 
been added. The descriptions below regarding test theory and the various cloud tests are 
applicable to GOES-10 data, as the difference in the last GOES-12 band’s wavelength 
requires some algorithm changes to the cloud mask that degrade accuracy of the cloud 
mask, principally in terminator regions (AFRL, 2001.) 

1, Remote Sensing Theory 

An understanding of basic satellite remote sensing theory is necessary before 
discussing the cloud tests, which rely upon various results from reflectance channels and 
emittance channels. Table 1.1 describes the satellite wavelengths used by the cloud mask 
tests, and the associated satellite system and resulting calibrated quantity. 
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Table 1.1. Sensor ehannel wavelengths used by the AER eloud mask algorithms. 
Wavelength Designations: VIS-visible, NIR-near infrared, SWIR-shortwave infrared, 
MWIR-mid-wave infrared, WV-water vapor, LWIR-long wave infrared. EBBT: 


Equivalent Blackbody Brightness Temperature. 


Tom AERL-VS-TR-2001-1549, p. 3. 


Wavelength 

(pm) 

Designation 

GOES 

Imager 

AVHRR 

MODIS 

Calibrated 

Quantity 

0.65 

VIS 

X 

X 

X 

Reflectance 

0.85 

NIR 


X 

X 

Reflectance 

1.38 

SWIR 



X 

Reflectance 

1.60 

SWIR 


X 

X 

Reflectance 

3.70-3.90 

MWIR 

X 

X 

X 

EBBT 

6.70-7.00 

WV 

X 


X 

EBBT 

8.55 

LWIR 



X 

EBBT 

10.8 

LWIR 

X 

X 

X 

EBBT 

12.0 

EWIR 

X 

X 

X 

EBBT 


a. Reflectance Channels 

Refleetanee properties for different surfaees vary with wavelength in the 
visible and near-infrared (NIR) bands. Water surfaces are poor reflectors in the NIR 
portion of the electromagnetic spectrum, whereas land surfaces are poor reflectors in the 
NIR except for at approximately 0.86 pm. Clouds, on the other hand are generally good 
reflectors in the NIR, except for ice clouds, which do not reflect well at 1.6 pm. Ice and 
snow exhibit peak reflectance at 0.65 pm and 0.85 pm. Taken together, these 
observations allow for automated discrimination between cloudy areas and clear areas. 

b. Emittance Channels 

In addition to using the NIR portion of the electromagnetic spectrum for 
cloud tests, the cloud algorithm also uses infrared (IR) channels as shown in Table 1. At 
these wavelengths, energy measured by the satellite sensor is mostly due to the emittance 
and temperature of the surface or atmospheric constituent (e.g., clouds), the transmittance 
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of the atmosphere. Refleetanee of eloud and background surfaces is negligible at these 
wavelengths. 

It is important to note that the 12.0 pm channel is more strongly affected 
by atmospheric water vapor absorption than the 10.8 pm channel, and as a result the 
brightness temperatures obtained with this channel are slightly less than those obtained 
from the 10.8 pm channel. Also, we note that in the 8.55 to 10.8 pm bands, absorption 
by atmospheric water vapor is roughly the same but absorption due to ice particles 
increases with wavelength between 10.8 and 12.0 pm. As we will see, these are all 
important factors used by the cloud mask algorithm in determining cloud presence and 
phase. 

c. Midwave Infrared Channel 

The 3.9 pm channel is unique because it provides both reflected solar 
energy and emitted energy. This characteristic, along with the optical properties of 
clouds at this wavelength, is very useful for detecting clouds that have weak or 
ambiguous spectral signatures at the other wavelengths used by the cloud detection 
algorithm. 

d. Water Vapor Channel 

We know that strong water vapor absorption occurs in the atmosphere at 
6.7 pm. Low-level features are masked out by this, resulting in resolution of only mid 
and high-level features. When cirrus optical depth is very small, transmission of 
background radiation through the cloud can reduce the cloud-background contrast to a 
level that is undetectable by automated cloud algorithms at visible, mid-wave IR, and 
long-wave IR wavelengths. This is especially true over bright and warm backgrounds 
such as deserts. Using data from the “opaque” 6.7 pm band, it is possible to develop a 
technique (see cloud test section below) to improve detection of thin cirrus over 
background conditions that stress traditional LWIR and visible techniques. 

2, Cloud Tests 

The cloud tests and background tests make use of the aforementioned satellite 
theory to determine if a given pixel is clear or cloudy, and the phase (ice or water) of 
cloudy pixels. Each test is performed by the cloud mask algorithm before final 
determination of cloud presence and phase is determined. Results from all tests are 
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considered jointly when assigning a elear or eloudy elassifieation and eloud phase (if 
applieable) to a given pixel. The eloud detection tests are described below: 

a. Thermally Distinct Cloud Test 

Most clouds that are not very low in the atmosphere exhibit a significantly 
colder brightness temperature than their terrestrial baekground. The thermally distinet 
eloud test is a simple single LWIR ehannel threshold test that is designed to diseriminate 
between the thermal signatures of obvious mid- and high-level elouds from the terrestrial 
baekground’s thermal signature. 

The test is not useful for highly transmissive cirrus (which can appear 
warmer than their true temperature), marine stratus and other low elouds, partially eloud- 
filled fields of view, and high-latitude snow/iee eovered backgrounds. Each of these 
cases is addressed by other eloud tests. 

This eloud test requires an estimate of the skin temperature over the region 
of interest so that eorresponding elear-sky brightness temperatures ean be eomputed. 
This skin temperature estimate is provided by a mesoseale model sueh as the AFWA 
MM5. As appropriate, a land-iee temperature model or a sea-surfaee temperature model 
is used to provide temperature data for iee and sea surfaees. 

The eloud test requires that the LWIR brightness temperature at 10.8 or 
12.0 pm be less than the predieted elear-seene temperature by an amount greater than the 
uneertainty in the clear seene estimate. The magnitude of the threshold is determined by 
aeeumulating statisties on the measured differenee between satellite EBBTs and the skin 
temperature of pixels determined to be eloud-free by previous runs of the cloud 
algorithm. Statistics are kept in a rotating 10-day database aeeording to satellite overpass 
time, loeation, and baekground geography type. 

b. Single Channel Reflectance Test 

This test diseriminates relatively high eloud refleetanee from a predieted 
baekground value in a similar manner to the thermal threshold test deseribed above. 
Clouds typieally have a higher refleetanee than their terrestrial baekgrounds at NIR 
wavelengths. To avoid diffieulties assoeiated with high refieetanee values over some 
land surfaees at 0.86 pm and water surfaees under eertain eonditions at 0.65 pm, the 


10 



scheme eonsiders only the results from the 0.86 |j,m band over water and the 0.65 [am 
band over land. In this manner, cloud contrast with the background surface is 
maximized. Separate thresholds are maintained for land baekgrounds and water 
baekgrounds. Limitations for this test inelude transmissive cirrus, high-albedo 
eeosystems, sun glint over water surfaees, desert, and snow/ice eovered baekground 
eonditions. 

c. Near-IR/Visible Reflectance Test 

This test is a ratio test that eompares the relative magnitudes of 0.86 pm 
(NIR) and 0.65 pm (visible) refleetanee values. Optically thick clouds generally exhibit a 
very similar speetral signature between the visible and NIR bands. For land and water 
surfaees, however, the speetral signatures between these two bands ean vary signifieantly. 
The ratio of NIR to VIS refleetanee can tell us a great deal about eloud eover. 

A ratio of approximately 1.0 indicates the likely presenee of eloud eover 
beeause surfaee features are effeetively masked. A ratio that differs signifieantly from 
1.0 indieates either vegetated land (in the ease of a ratio greater than 1.0) or water (in the 
ease of a ratio less than 1.0.) The thresholds between elouds/no elouds are determined 
empirieally, and are adjusted as needed when eonditions sueh as high humidity and high 
aerosol/haze eoneentration exist. The results of this test are only useful in the absenee of 
sun glint, desert baekground, and snow/iee eovered baekground eonditions. Additionally, 
mixed geography types sueh as eoastlines prevent use of this test. 

d. Low Cloud and Fog Test 

Low elouds and fog are typically virtually indistinct from their 
baekground surfaees in the LWIR, but show high contrast to their backgrounds at visible 
and NIR wavelengths during daytime. Low elouds and fog also have recognizable 
speetral signatures in the MWIR when eompared to LWIR brightness temperatures. The 
low eloud and fog test is designed to supplement the visible and NIR test during daytime, 
and addresses the inability of the thermal eloud test to diseriminate thermally indistinet 
elouds at night. 

The magnitude of the MWIR-LWIR speetral differenees for small water 
droplet elouds varies depending on time. During daytime, a eloud’s refleetivity in the 

3.9 pm band is greater than the LWIR brightness temperature for the eloud. At night, the 
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emissivity of a cloud measured in the 3.9 |im band is less than the LWIR-indieated 
brightness temperature. The sign of the signature ehanges from day (positive differenee) 
to night (negative differenee.) 


This test maintains separate thresholds for desert, non-desert, and potential 
sun glint backgrounds due to the high refleetance within the 3.9 pm band of those 
features. Desert surfaees have a lower 3.9 pm emissivity at night than other terrestrial 
surfaces, so separate nighttime thresholds are maintained for desert areas. 

3, Cirrus-Specific Cloud Tests 

As the eloud deteetion tests have shown, eirrus (particularly thin cirrus) are often 
difficult to detect. Optically thick cirrus clouds are easily identified via temperature 
eontrast with the baekground. However, cirrus clouds with nonzero transmittanee are 
common, resulting in LWIR brightness temperatures that can be mueh greater than the 
actual cloud temperature. This ean result in eirrus being erroneously elassified as mid¬ 
level (water) elouds by automated classifieation sehemes. The sehemes deseribed below 
help to improve deteetion of cirrus and ensure eorreet phase categorization. The AFRL 
(2001) technieal report was used as a referenee in preparing this seetion, and eontains 
further details regarding eaeh of the tests deseribed below. 

a. LWIR Cirrus Test 

This test examines LWIR brightness temperature differences between 
8.55, 10.8, and 12.0 pm. As discussed earlier, iee particle absorption inereases with 
wavelength aeross this speetral region and water vapor attenuation is slightly stronger at 
the 12 pm wavelength. These faetors result in eirrus brightness temperatures that 
deerease with wavelength. 

The eloud deteetion threshold used for this test is defined as a function of 
atmospheric water vapor and path length through the atmosphere due to the faet that, 
given the faetors described above, it is possible to obtain erroneous results due to 
atmospherie moisture in a elear scene. Additionally, over snow or iee-eovered 
baekgrounds, the test assumes that 10.8 pm brightness temperatures measured from the 
eirrus clouds are colder than the clear-seene terrestrial baekground by an amount greater 
than a defined threshold value. 
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b. NIR Cirrus Test 

Optically thin cirrus can be detected using the Moderate Resolution 
Imaging Speetroradiometer (MODIS) in the 1.38 pm band. Here, water vapor absorption 
masks out refleeted solar radiation from beneath elouds so that the eontrast between 
cirrus clouds and the baekground is suffieient to allow for thin eirrus deteetion. 

This test has diffieulty over areas with high terrain relief (snow-eovered 
mountains) and loeations with very low atmospheric water vapor and refleetive surfaees 
(polar areas.) 

c. MWIR-LWIR Cirrus Test 

The 12.0 pm ehannel is typieally used for eomparison with the 3.9 pm 
ehannel in this test beeause of ice partieles’ inereased absorption at this wavelength 
eompared to the 10.8 pm wavelength. In practiee, however, there is little differenee in 
the observed temperatures between the 10.8 pm and 12.0 pm wavelengths. As a result, 
the eloud signature obtained by differeneing the 3.9-12.0 pm bands is only slightly less 
than the signature obtained by differeneing the 3.9-10.8 pm bands is only slightly greater 
for the 3.9-12.0 pm band. 

Use of the 12.0 pm ehannel, however, ean result in errors that lead to false 
eloud detection when there is a clear seene with high humidity levels near the surface. In 
this situation, strong water vapor absorption near the surfaee eauses the 12.0 pm 
equivalent blaekbody temperature (EBBT) to appear mueh lower than the true 
temperature. To remedy this problem, the 10.8 pm data are used rather than the 12.0 pm 
data when integrated water vapor values in a eolumn are greater than a defined threshold 
value. 

d. Water Vapor Channel Cirrus Test 

The 6.7 pm wavelength measures upwelling thermal radiation emitted by 
water vapor in the atmosphere. The magnitude of the 6.7 pm brightness temperature 
depends on the amount of water vapor present along the atmospherie path, but will be 
less than temperatures eomputed for the same seene at LWIR wavelengths in all but the 
driest of atmospheres. Any eloud in the field of view will reduee the differenee between 
the 6.7 and 10.8 pm brightness temperatures sinee 10.8 pm brightness temperature will 
be representative of the ambient temperature at the eloud’s level. 

13 



The water vapor ehannel eirrus test first determines known elear-seene 
areas, evaluates the 6.7-10.8 pm temperature differenees for these areas, and then 
searehes the seene for loeations where the differenee falls below a eutoff threshold based 
on the elear-seene differenee. This test is advantageous beeause it does not require us to 
know surfaee eharacteristies, whieh in some eases may be difficult to obtain. 

4. Cloud Phase Classification 

In order to obtain cloud optical and microphysical properties, it is necessary for 
the cloud mask algorithm to discriminate between ice cloud and water cloud. The cloud 
tests above use spectral signatures that can represent either ice cloud or water cloud, so 
phase discrimination is still necessary. 

The phase test uses spectral signatures representative of ice cloud to classify the 
phase of pixels identified as cloudy by the aforementioned tests. If a pixel fails to meet 
ice cloud spectral signature requirements, it is classified as a water cloud pixel. The 
phase determination is made by applying several tests, as described below. 

a. LWIR Brightness Temperature Test 

Pixels are tested to see if their brightness temperature value allows them to 
be unequivocally classified as ice cloud. Any pixel with a brightness temperature less 
than 233K is classified as ice cloud. 

b. MWIR-LWIR Test 

A positive MWIR-LWIR difference (at night) is unique to transmissive 
cirrus for reasons described above. 

c. NIR Test 

During daytime, if the NIR thin cirrus test detects cloud and the other 
reflectance tests do not, then the pixel is classified as ice cloud. 

d. Brightness Temperature Difference Test 

Due to the nearly uniform increase in particle absorption with wavelength 
in the LWIR, coupled with the preferential increase between 10.8 and 12.0 pm for water 
droplets, an 8.55-10.8 pm BTD greater than a 10.8-12.0 pm BTD difference is indicative 
of ice cloud, while an 8.55-10.8 pm BTD less than a 10.8-12.0 pm BTD is indicative of 
water cloud. 


14 



C. AIR FORCE WEATHER AGENCY MESOSCALE MODEL 5 

The Mesoscale Model 5 (MM5) was developed by the National Center for 
Atmospherie Research (NCAR) and the Pennsylvania State University. The MM5 was 
declared operational for Air Force use in 1997 by the Air Force Weather Agency 
(AFWA), and is tailored as necessary by AFWA to perform as well as possible in every 
theatre in which it is run. Various physics packages determine how the model treats 
various atmospheric processes such as convection, cloud moisture phase changes, etc. 
Certain parameters contained in the physics packages can be adjusted by modelers at 
AFWA to ensure the most favorable configurations are used. 

The MM5 implementation used in this thesis is the 15km AFWA MM5. This 
model contains 42 sigma levels, with highest vertical resolution near the surface and 
decreasing vertical resolution as we travel toward the top of the modeled atmosphere at 
50mb. This is a critical factor when considering cirrus forecasting within the model. 
Since cirrus are often only a few meters thick, their location and ice water concentration 
may not be accurately represented by the model. Various papers have studied the 
problem of forecasting cirrus using a mesoscale model. Starr and Wylie (1990) found 
that vertical resolution of 0.5 km or better is probably required to adequately resolve 
dynamical forcing required for the development of a physically-based model 
parameterization of cirrus. 

The problem of forecasting cirrus within a numerical model is an extremely 
difficult and complex task. Forecasting of cirrus goes well beyond simply knowing 
where cirrus are initially located and advecting the clouds properly. Not only must the 
model account for the movement and development of cirrus, it must also account for 
development related to microphysical processes that occur many orders of magnitude 
below the model’s resolution. These cloud microphysical processes are represented 
within a forecast model by any one of several parameterizations. The AFWA MM5 uses 
a modified Reisner mixed-phase scheme. Due to the tendency to overforecast ice, the 
scheme as implemented in the AFWA MM5 does not account for riming processes or 
graupel processes when computing ice water mixing ratio. 
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Norquist and d’Entremont. (2003) have used explicit ice water content to 
represent cirrus in MM5 output. The ice water mixing ratio field output by the model is 
the result of many calculations for each model gridpoint, and is dependent on several 
factors. Equation 1 shows the ice water mixing ratio formula as given by Reisner: 


dp * qi 


= -ADV{P*q,) + DIV{P ^q.) + D{q,) + 


D*(P +P +P +P +P -P -P -P -P -P 
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It is convenient to view handling of ice water mixing ratio in the model in terms of 
sources and sinks (Reisner, 1998). Table 2.1 categorizes the production terms in the 
equation above into sources and sinks. Items listed in Reisner’s equation that are not 
included in the AEWA MM5 are Picng, Pi.iacw, and Praci- 


Table 2.1. Terms comprising the Reisner mixed-phase ice microphysics scheme. 
Shaded areas indicate scheme terms that do not apply to the AEWA MM5 
implementation of the scheme. __ 


Source Terms 

Sink Terms 

Pidsn: Initiation of nucleation 

Psaci: Snow generation due to collection of 

ice by snow 

Pifzc: Heterogeneous and homogeneous 

freezing of cloud drops 

Picns: Total conversion of cloud ice to snow 

Pispi: Increase in mass of cloud ice 

associated with ice multiplication 

Pimit: Cloud water generation by the 

melting of ice 

Pidep: Depositional growth of cloud ice 

Picng: Conversion of cloud ice to graupel 

Pi.iacw: Contribution of riming to ice 

production 

Praci: Graupel generation by collision of ice 

and cloud water 


Reisner expresses initiation of ice nucleation as a function of the number of ice 
crystals forecast at a particular gridpoint versus the number of ice crystals already present 
at that gridpoint. The initial number of crystals at a gridpoint is empirically determined 
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using Fletcher’s curve (Fletcher, 1962) and temperature data taken from the model. The 
scheme is known to overpredict ice crystal concentrations for very low temperatures, so 
the lowest temperature is limited to 246 K. If the number of ice crystals predicted by the 
Fletcher curve is greater than the number of ice crystals already present at a given 
gridpoint, the additional vapor at that gridpoint is used to generate ice crystals until the 
predicted number and actual number match, at a rate determined by the first source term 
listed in the table above. With a known number of ice crystals, the model can now 
calculate the rate of vapor depositional growth of cloud ice, which depends on the extent 
of supersaturation, mean crystal size or mass, and the number of ice crystals. Vapor 
transfer rates for riming and graupel processes are neglected by the AFWA MM5. 
Additionally, the MM5 normalizes vapor deposition rates so that over depletion of vapor 
and subsaturation (associated with latent heat release) do not take place. 

The model computes the results for each process identified in the table above 
based upon formulas (reference NCAR Technical Note TN-398+STR for details), and 
uses those results as part of its process to determine results for the ice water mixing ratio 
field. There are many uncertainties involved in the computation of ice water values for a 
given location. For instance, determination of the number of particles at a given location, 
based on statistically-determined particle concentration, is a source of error for 
heterogeneous nucleation. Temperature thresholds for homogeneous nucleation are 
uncertain, and this temperature uncertainty can contribute to errors in the forecast number 
of particles at a gridpoint. Minus 40 C is the value used by the MM5 to represent the 
onset of homogeneous nucleation; however homogeneous nucleation may occur at 
slightly warmer temperatures around minus 38 C in the real atmosphere. 
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III. METHODS AND PROCEDURES 


A. ANALYSIS METHOD 

Cases selected for this study are typical for day-to-day forecasting operations. 
Cases were selected by observing infrared satellite loops for periods of significant cirrus 
field development and/or evolution. The focus for each case is the 18-30 hour forecast 
period, so cases were selected so that the time of the initial image coincided with the 18- 
hour forecast time of the MM5 run that was valid at the time the case began. 

The MM5 provides forecast output for each 3-hour point within the period of 
interest, so satellite and model data were also compiled at times corresponding to the 
MM5’s 18, 21, 24, 27, and 30 hour forecast periods. The Naval Postgraduate School 
Meteorology department typically receives GOES satellite data via FTP from the Naval 
Research Laboratory, Monterey and 15km MM5 pressure and sigma level data are 
routinely received from the Air Force Weather Agency for select forecast hours from 
each model run. 

Once all satellite images for each case were saved, the ABF cloud mask 
processing code was applied to the images and merged with the image data. Proprietary 
scripts were used to place these data into a form useful by other routines that arranged the 
data into a format that could be studied along with model output within GEMPAK/GARP 
(General Meteorology Package / General Analysis and Rendering Program) software. 

B. CASE STUDY MODEL VERIFICATION METHOD 

In order to gain an understanding of how the MM5 forecast predicted the basic 
atmospheric structure from the case study start time through the 30-hour forecast, the 
observed Global Forecast System (GFS) analyses were used to verify the intensity, 
placement, and movement of significant features. The GFS analyses were chosen to 
verify the MM5 due to their increased availability versus subsequent MM5 model 
analyses which are only available every 12 hours. This allowed for direct verification of 
the MM5 more frequently during the study period for each case. Height and wind fields 
at various levels were considered in the verification and only areas where significant 
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deviation between the GFS and MM5 oeeurred are shown in the figures assoeiated with 
eaeh ease. 

To assess the model dynamies, moisture, and eloud distributions, the MM5 
foreeasts were eompared to the satellite-derived eloud mask. This was done using GARP 
to overlay relevant fields. Position and phase errors in the model foreeasts—partieularly 
with the first ease study—prevented direet eomparison of the model iee field to observed 
eirrus elouds, but the eharaeteristies of the model elouds and assoeiated dynamie foreings 
were still useful to help understand how the model generated the elouds. 
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IV. CASE STUDY I: UPPER-LEVEL CLOSED LOW 


A. CASE DESCRIPTION 

This case focused on a elosed upper-level low moving toward the western United 
States on 21 October 2005. The low had no jet support, and moved very slowly to the 
east during the study period. Figures 4.1.a through 4.1.e show the evolution of the cirrus 
field from 21/OOZ through 21/12Z. The cirrus field at the beginning of the ease study 
was mostly loeated over the Paeifie, and was the result of eonveetion over the tropies the 
day before the ease study began. The cirrus field subsequently entered an area of 
inereasing difiuenee between an upper low to its north and a ridge to its east that aeted to 
spread out the eirrus at the northern extent of the plume. Additionally, large-scale 
upward motion in the difiuent region resulted in the development of new cirrus. Some of 
the eirrus clouds were pulled downstream, to the east, around flow related to the ridge. 
Other eirrus elouds were drawn around the upper low by the end of the period. 

B, MM5 FORECAST VERIFICATION 

The GFS analyses from 21 Oetober at OOOOZ, 0600Z, and 1200Z were used to 
verify the eorresponding 18, 24, and 30-hour foreeasts from the 20 Oetober 0600Z MM5 
run. Verification was performed for wind and height fields at 300, 400, and 500mb to 
determine how well the model replieated the cirrus-level transport winds. 

Generally, foreeast errors in the wind speed field were minimal but signifieant 
errors were present in the height field by the 18 hour foreeast period. The errors present 
in the height field at the 18 hour foreeast period persisted through the end of the ease 
study. Figure 4.2.a shows the 18-hour MM5 height foreeast verifieation at 300mb. The 
height of the low is deeper in the GFS analysis than for the MM5 forecast, and the ridge 
over the interior of Mexieo is represented more strongly in the GFS than the MM5. Most 
signifieantly, however, the MM5 misses the sharpness of the ridge axis aeross Mexieo 
and the desert southwest. The result of these height errors is to produee flow that is more 
zonal in the MM5 output than is observed. Consequently the aetual eloud field (shown in 
Figures 4.1.a through 4.1.e) appears to agree better with the GFS analyzed height field 
with its more north-south orientation. In addition, direetly ahead of the low, the GFS 
indieates mueh stronger difiuenee taking plaee than is oeeurring in the MM5. This 
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analyzed difference in the GFS is consistent with the fanning out of the cirrus field shown 
in Figures 4.1.C through 4.1.e. While the MM5 shows some difluence, it tends to be less 
and is located in a different area than the observed cirrus. Therefore, the MM5 cirrus 
forecasts associated with this process are characterized by significant placement error. 

The MM5 wind field (not shown) indicated a jet feature around the base of the 
low that extends through southern Arizona, central New Mexico, and across northern 
Texas. This band of higher winds is located just to the north of the modeled cirrus 
coverage, and is typical of cirrus forced by an anticyclonically curved jet. 

Figure 4.2.b shows the 24-hour MM5 forecast compared to the GFS analysis valid 
at 06Z on 21 Oct 05. The height pattern remains similar to what is shown in Figure 4.2.a, 
with the GFS showing a much more marked ridge across the desert southwest than the 
MM5. The height discrepancies are more evident at mid-levels for this forecast period, 
which indicates the MM5 is beginning to deviate more from observed conditions. The 
MM5 continues to underforecast the depth of the low near the west coast and intensity of 
the ridging east of the low. The differences in extent of the ridge axis are particularly 
evident across central Arizona into southern Nevada, which is where the observed cirrus 
shown in Figure 4.1.C is beginning to spread out. 

As Figure 4.2.c shows, the GFS continues to analyze a deeper low at 300mb at 
21/12Z than the 30-hour forecast from the MM5. This trend was found at lower levels as 
well. The GFS also continues to analyze a sharper ridge across the desert southwest and 
northern Mexico than forecast by the MM5. The GFS and MM5 placements of the low 
center between 300mb and 500mb (not shown) agree very well, as has been the case 
throughout the entire period of interest. This indicates that the MM5 is moving the low 
appropriately despite weak synoptic-scale support for its movement. Given the GFS 
upper air height pattern, the analyzed GFS winds at 200mb, 300mb, and 500mb are from 
a more southerly direction ahead of the low as noted earlier, and the height pattern overall 
remains more conducive to stronger divergence (and greater cloud coverage) ahead of the 
low than the forecast by the MM5 indicates. 
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C. ANALYSIS 

The upward motion and ice water mixing ratio were examined to determine the 
relationship between the dynamics and the model-produced cirrus clouds. Case I is 
characterized by strong upward vertical motion aloft in the difluent area ahead of the 
closed low, and strong upward motions streaming eastward across northern Mexico, as 
shown in Figures 4.3.a through 4.3.e. The main band of ice water mixing ratio values 
began at approximately 1725 sigma and continued to about 3190 sigma (approximately 
170mb to 320mb), although there were occasional small areas of ice water mixing ratio 
values above 1725 sigma. The ice water field gradually diminished below 3190 sigma. 

The 2880 sigma level, approximately 280mb, was chosen as the representative 
sigma level for the modeled cirrus field for this case. Normally cirrus are expected in 
areas of fairly weak upward motion, and indeed there is a strong correlation between 
upward motion and the presence of ice water in the model at 2880 sigma as indicated in 
Figure 4.3.a through 4.3.e. (Note: vertical motions are taken from the 300mb level but 
are representative of vertical velocities within the main cirrus band at all levels where it is 
found.) However, there is also a strong correlation between the 2880 sigma cirrus field 
and moisture divergence, as shown in Figures 4.4.a through 4.4.e. The correspondence 
between moisture divergence and presence of ice water mixing ratio values at the 2880 
sigma level shows that the model ice must be continually created to offset the moisture 
divergence. The creation process is the upward motion from below of nearly saturated 
air (with respect to ice) that freezes into ice as it cools. Critical to maintaining cirrus of 
this type in the model is the supply of nearly saturated air in the slowly rising plume. If 
this diminishes, the model clouds may dissipate too soon. 

Interestingly, subjective analysis showed that the amount of ice coverage at the 
upper sigma levels in and near the top of the band diminished with time, particularly after 
the 06-hour forecast, while the amount of ice coverage at lower levels increased. The 
system does become less organized with time, and this is likely a reflection of 
diminishing support for cirrus formation as upward motions and, consequently, support 
for large-scale supersaturation subside. The decrease in organization of model ice also is 
evident in the observed cirrus fields (Figures 4.1.a through 4.1.e.) 
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In general, the iee water tends to coineide with regions of relatively strong upward 
motion (shaded areas in Figure 4.3.a through 4.3.e.) The broad area of upward motion 
ahead of the low over northern Baja produees no ice water in the model. The lack of ice 
water in this region is presumably due to the fact that the moisture source for the ice 
water does not flow into this region. 

Additionally, Figures 4.5.a through 4.5.e indicate that the upward vertical velocity 
fields lie within areas of strong frontolysis. Again, this was an unexpected result. As 
time passes, however, the cirrus field dissipates due to weakening of the low and ongoing 
frontolytical forcing. It appears that this frontolytical forcing provides an inverse thermal 
circulation, which results in slight warming of air within the cirrus field and therefore 
causes a decrease in relative humidity with respect to ice over time. 

D, RESULTS 

To verily the model’s cirrus forecast, the ice water field composite for all levels 
where ice water was present in the model was plotted versus the corresponding cloud 
mask results from 21/OOZ through 21/12Z as shown in Figures 4.6.a through 4.6.C. It is 
evident that the ice water field is dissipating with time, which agrees well with the 
frontolytical forcing evident in Figures 4.5.a through 4.5.e. Significant phase errors were 
present from the 18 hour forecast through the 30 hour forecast. As described earlier, this 
was a result of the MM5 under-forecasting the ridge sharpness, and as a result under¬ 
forecasting the location of cirrus as well as the strength of divergent flow across Arizona 
which led to significant cirrus coverage there. Finally, overall cloud coverage was less 
than observed via the cloud mask results, particularly later in the forecast period when the 
MM5 began to forecast diminished frontogenesis and upward vertical motions along the 
main cirrus band over northern Mexico and the southern United States. 

The ice field shown in Figures 4.7.a through 4.7.C was examined to determine 
whether the underforecast of model cirrus might be due to a lack of moist air to feed the 
clouds. Again, our ice water field was becoming more diffuse with time, but the satellite 
images in Figure 4.1.a through 4.1.e show that the observed cirrus pattern remained 
relatively intact with time. Importantly, the 90% relative humidity over ice (shaded 
green) appears to fit the observed cirrus pattern better than the ice water mixing ratio 

composite when model phasing error is taken into account. This suggests that the 
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moisture supply was sufficient but the dynamic forcing was not sufficiently strong to 
produce enough cirrus in the model. The model may require too much vertical motion to 
initiate and maintain ice clouds, which contributes to its underforecast of cirrus that is 
observed in the real atmosphere accompanying far lower vertical velocities. 
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Figure 4.1.a. GOES-W image valid at 21 Get 05/00Z. 



Figure 4.1.b. GOES-W image valid at 21 Get 05/03Z. 
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Figure 4.1.c. GOES-W image valid at 21 Oct 05/06Z. 


Figure 4.1.d. GOES-W image valid at 21 Oct 05/09Z. 
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Figure 4.1.e. GOES-W image for valid at 21 Oct 05/12Z. 



Figure 4.2.a. Verification of the 18-hour MM5 300mb height forecast 
(green contours) against the corresponding GFS analysis (blue contours) 
for OOZ on 21 Oct 05. 
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Figure 4.2.b. Verification of the 24-hour MM5 300mb height forecast 
(green contours) against the corresponding GFS analysis (blue contours) 
for 06Z on 21 Oct 05. 



Figure 4.2.c. Verification of the 30-hour MM5 300mb height forecast 
(green contours) against the corresponding GFS analysis (blue contours) 
for 12Z on 21 Oct 05. 
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Figure 4.3.a. 300mb vertical motion and 2880 sigma ice water mixing 
ratio for 21/OOZ. Ice water mixing ratio is in grams per kilogram x 10' , 
and upward vertical motion (represented by warm shades) is in m s'*. 



Figure 4.3.b. 300mb vertical motion and 2880 sigma ice water mixing 
ratio for 21/03Z. Ice water mixing ratio is in grams per kilogram x 10'^, 
and upward vertical motion (represented by warm shades) is in m s'*. 
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Figure 4.3.c. 300mb vertical motion and 2880 sigma ice water mixing 
ratio for 21/06Z. Ice water mixing ratio is in grams per kilogram x 10' , 
and upward vertical motion (represented by warm shades) is in m s'*. 



Figure 4.3.d. 300mb vertical motion and 2880 sigma ice water mixing 
ratio for 21/09Z. Ice water mixing ratio is in grams per kilogram x 10'^, 
and upward vertical motion (represented by warm shades) is in m s'*. 
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Figure 4.3.e. 300mb vertical motion and 2880 sigma ice water mixing 
ratio for 21/12Z. Ice water mixing ratio is in grams per kilogram x 10'^, 
and upward vertical motion (represented by warm shades) is in m s'*. 



Figure 4.4.a. 2880 sigma ice water mixing ratio and moisture 
divergence for 21/OOZ. Ice water mixing ratio is in grams per kilogram 
and areas of moisture divergence are shaded. 


32 













Figure 4.4.b. 2880 sigma ice water mixing ratio and moisture 
divergence for 21/03Z. Ice water mixing ratio is in grams per kilogram 
and areas of moisture divergence are shaded. 



Figure 4.4.c. 2880 sigma ice water mixing ratio and moisture 
divergence for 21/06Z. Ice water mixing ratio is in grams per kilogram 
and areas of moisture divergence are shaded. 
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Figure 4.4.d. 2880 sigma ice water mixing ratio and moisture 
divergence for 21/09Z. Ice water mixing ratio is in grams per kilogram 
and areas of moisture divergence are shaded. 



Figure 4.4.e. 2880 sigma ice water mixing ratio and moisture 
divergence for 21/12Z. Ice water mixing ratio is in grams per kilogram 
and areas of moisture divergence are shaded. 
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Figure 4.5.a. 2880 sigma ice water mixing ratio and frontogenesis for 
21/OOZ. Ice water mixing ratio is in grams per kilogram, and 
frontogenesis is represented by shaded areas. 



Figure 4.5.b. 2880 sigma ice water mixing ratio and frontogenesis for 
21/03Z. Ice water mixing ratio is in grams per kilogram, and 
frontogenesis is represented by shaded areas. 
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Figure 4.5.c. 2880 sigma ice water mixing ratio and frontogenesis for 
21/06Z. Ice water mixing ratio is in grams per kilogram, and 
frontogenesis is represented by shaded areas. 



Figure 4.5.d. 2880 sigma ice water mixing ratio and frontogenesis for 
21/09Z. Ice water mixing ratio is in grams per kilogram, and 
frontogenesis is represented by shaded areas. 
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Figure 4.5.e. 2880 sigma ice water mixing ratio and frontogenesis for 
21/12Z. Ice water mixing ratio is in grams per kilogram, and 
frontogenesis is represented by shaded areas. 



Figure 4.6.a. Composite ice water mixing ratio field vs cloud mask 
result for 21/OOZ. Ice water mixing ratio returns are contoured, and ice 
cloud results from the cloud mask are in blue. 
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Figure 4.6.b. Composite iee water mixing ratio field vs eloud mask 
result for 21/06Z. Iee water mixing ratio returns are contoured, and ice 
cloud results from the cloud mask are in blue. 



Figure 4.6.C. Composite ice water mixing ratio field vs cloud mask 
result for 21/12Z. Ice water mixing ratio returns are contoured, and ice 
cloud results from the cloud mask are in blue. 
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Figure 4.7.a. Relative humidity (with respect to ice) vs composite ice 
water mixing ratio for 21/OOZ. Green shades indicate relative humidity 
of 90% and greater. 


Figure 4.7.b. Relative humidity with (respect to ice) vs composite ice 
water mixing ratio for 21/06Z. Green shades indicate relative humidity 
of 90% and greater. 
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Figure 4.7.c. Relative humidity (with respeet to iee) vs composite ice 
water mixing ratio for 21/12Z. Green shades indicate relative humidity 
of 90% and greater. 
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V. CASE STUDY II: UPPER-LEVEL RIDGE OVER ROCKIES 


A. CASE DESCRIPTION 

Case II involved eirrus assoeiated with an upper-level trough/ridge pattern over 
the Roeky Mountain States from 26 Oct 05/00Z through 26 Oct 05/12Z. Figures 5.1.a 
through 5.1.e show the GOES infrared satellite images every three hours during the case 
study valid time. At 26/OOZ, the main areas of interest were an elongated band of cirrus, 
which diminishes with time, located across the northern and central plains states and a 
more loosely organized area of ice cloud located over Utah. The ice clouds over Utah 
were the result of lift around a decaying upper low located in southeastern Nevada just 
prior to case study start time, whereas the cirrus clouds over the plains states were the 
result of advection across a ridge axis of moisture originally lifted by the low. Due to the 
configuration of the data collection system at the time the data were obtained for this 
case, 21-hour sigma data are not available for this case study. 

B, MM5 FORECAST VERIFICATION 

The 300mb height field for the case study valid time is shown in Figures 5.2.a 
through 5.2.C. The placement and intensity of the trof at 300mb is very similar between 
the GFS analysis and the MM5 forecast at all times, but the analyzed versus forecast 
ridge positions and intensities show minor differences. GFS-analyzed winds (not shown) 
were 10 to 15 knots stronger downwind of the ridge axis between the U.S.-Canada border 
and the central plains. This difference in wind speed was present at all analyzed levels 
between 200mb and 500mb. 

One important difference between the GFS analyses and the MM5 forecasts is 
that the GFS appears to indicate a shortwave trof progressing through the pattern to the 
east of our area of interest, with a corresponding change in alignment of the height 
contours and resulting geostrophic winds. This results in a slightly more confluent area 
aloft over the northern and central plains, which may play a role in how the cirrus field in 
that area evolves during the forecast period. 
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C. ANALYSIS 

Case Study I indicated that the presence of ice water in the model was tied 
strongly to upward vertical motion. This is valid in this case study as well, even in the 
area downstream of the ridge where production should not occur, as evident in 
comparisons of the upward motions shown in Figures 5.3.a, 5.3.C, and 5.3.e to the ice 
water field for corresponding times shown in Figures 5.4.a through 5.4.c. The vertical 
motion fields shown in Figures 5.3.a through 5.3.e indicate very weak, large-scale 
upward motion associated with the area downwind of the trof over Nevada and 
immediately east of the ridge axis. These upward motions are not very well organized, as 
evident in the pockets of upward and downward motions over eastern Colorado, 
Wyoming, and Montana in Figure 5.3.a. The upward motion field east of the ridge axis 
becomes more organized by 06Z (as seen in Figure 5.3.c) as a result of the trough lifting 
northeastward in the MM5 (Figure 5.2.b.) Stronger, better organized upward motion is 
located near the upper low over Utah throughout the case study period. 

As mentioned earlier, this case provides two distinct areas of interest with regard 
to ice cloud coverage. For the first area, downstream of the ridge axis, the model’s ice 
returns were located at higher altitudes—generally above 4145 sigma. Peak ice water 
mixing ratio values downstream of the ridge axis occur at 1460 sigma as shown in 
Figures 5.5.a through 5.5.d. Ice water mixing ratio results indicate diminishing ice water 
coverage within this plume as time progresses, and agree well with observed cirrus 
location. The second distinct area of ice cloud coverage was closer to the upper-level 
low, upstream of the ridge axis. Here, the most extensive ice water coverage ahead of the 
low is concentrated between the 4145 and 6065 sigma levels, with some traces of ice 
remaining at lower levels. In this location ice water returns are not as consistent 
vertically as in the first case. Additionally, convective elements are present in the returns 
for this case as indicated by darker contoured areas resembling bullseyes in Figures 5.4.a 
through 5.4.C and 5.6.a through 5.6.C. The model overforecast ice water coverage in this 
area as evident in Figures 5.4.a through 5.4.c, particularly over Colorado, Wyoming, 
Idaho, and Utah. In these locations, abundant ice water was present in areas absent of 
observed clouds. 
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D, RESULTS 

With regard to basic forecast error, the second case posed somewhat less of a 
challenge for the MM5 than the first case. Figures 4.2.a through 4.2.c show significant 
phase error in the model versus observed conditions, whereas Figures 5.2.a through 5.2.c 
show the model handled evolution of the synoptic situation fairly well. Discrepancies did 
arise in the position of the downstream ice cloud in Figures 5.4.a through 5.4.c, indicating 
the presence of a slight timing error. Additionally, along this downstream plume, 
coverage of ice was more continuous in the model than indicated in the cloud mask. 
Near the low, significant differences arose in the placement of ice cloud. The model 
appeared to tie ice cloud very strongly to upward motions, so the presence of sufficient 
moisture and upward motion appeared to result in an overforecast of ice cloud in this 
location. Comparison of the ice water field to the cloud mask, shown in Figure 5.4.a, 
shows that ice cloud over northern Arizona was slightly underforecast, although the 
model did a slightly more respectable job of forecasting ice cloud over Utah at this time. 
Significant overforecasting of ice cloud occurred over western Colorado, southwestern 
Wyoming and eastern Idaho, corresponding to areas of upward motion as shown in 
Figure 5.3.a. Figure 5.4.b shows significant error across Wyoming, and also some timing 
and coverage errors associated with the downstream plume across northeastern Colorado. 
By 26/12Z, the model did a very respectable job depicting ice water coverage over 
Colorado, but appeared to underforecast ice over northern New Mexico. This appears to 
be due to insufficient moisture in the area, as mixing ratio values across the area were 
initially very low and increase only slightly by 26/12Z. 

The 90% relative humidity (with respect to ice) field shows nearly saturated air 
over Wyoming in Figure 5.6.a whereas the corresponding satellite image (shown in 
Figure 5.1.a) and cloud mask image (shown in Figure 5.4.a) show only a swath of high 
moisture over eastern Wyoming. The relative humidity field suggests that ample 
moisture to the northeast of the low, when coupled with vertical motion of sufficient 
strength, leads to the erroneous cloud cover over Colorado, Wyoming, and Utah 
mentioned above. The broad area of high relative humidity and loosely organized 
vertical motion lead to cirrus in areas that do not match the observed clouds in detail. 
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This represents a signifieant foreeast ehallenge to get the vertieal motion exaetly right so 
that eirrus are not over predicted in this broad area of high relative humidity. 




Figure 5.1.a. GOES-W image valid at 26 Oct 05/00Z. 


Figure 5.1.b. GOES-W image valid at 26 Oct 05/03Z. 
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Figure 5.I.C. GOES-W image valid at 26 Oct 05/06Z. 



Figure 5.1.d. GOES-W image valid at 26 Oct 05/09Z. 
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Figure 5.1.e. GOES-W image valid at 26 Get 05/12Z. 



Figure 5.2.a. Verifieation of the 18-hour MM5 300mb height foreeast 
(green contours) against the corresponding GFS analysis (blue contours) 
for OOZ on 26 Oct 05. 
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Figure 5.2.b. Verification of the 24-hour MM5 300mb height forecast 
(green contours) against the corresponding GFS analysis (blue contours) 
for 06Z on 26 Oct 05. 



Figure 5.2.c. Verification of the 30-hour MM5 300mb height forecast 
(green contours) against the corresponding GFS analysis (blue contours) 
for 12Z on 26 Oct 05. 


47 



Figure 5.3.a. 400mb vertical motion and 4145 sigma ice water mixing 
ratio for 26/OOZ. Upward vertical motion is represented by warm 
shades, and ice water mixing ratio is in blue contours. 



Figure 5.3.b. 400mb vertical motion at 26/03Z. Sigma level data were 
not available for this time; as a result, the ice water field is not available. 
Upward vertical motion is represented by warm shades. 
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Figure 5.3.c. 400mb vertical motion and 4145 sigma ice water mixing 
ratio for 26/06Z. Upward vertical motion is represented by warm 
shades, and ice water mixing ratio is in blue contours. 



Figure 5.3.d. 400mb vertical motion and 4145 sigma ice water mixing 
ratio for 26/09Z. Upward vertical motion is represented by warm 
shades, and ice water mixing ratio is in blue contours. 
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Figure 5.3.e. 400mb vertical motion and 4145 sigma ice water mixing 
ratio for 26/12Z. Upward vertical motion is represented by warm 
shades, and ice water mixing ratio is in blue contours. 



Figure 5.4.a. Composite ice water mixing ratio field vs cloud mask 
result for 26/OOZ. 
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Figure 5.4.b. Composite ice water mixing ratio field vs cloud mask 
result for 26/06Z. 



Figure 5.4.c. Composite ice water mixing ratio field vs cloud mask 
result for 26/12Z. 
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Figure 5.5.a. 1460 sigma ice water mixing ratio contours and 
corresponding IR image valid at 26/OOZ. Labels have been omitted from 
the ice water mixing ratio contours. 



Figure 5.5.b. 1460 sigma ice water mixing ratio contours and 

corresponding IR image valid at 26/06Z. Labels have been omitted from 
the ice water mixing ratio contours. 
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Figure 5.5.C. 1460 sigma ice water mixing ratio contours and 

corresponding IR image valid at 26/09Z. Labels have been omitted from 
the ice water mixing ratio contours. 



Figure 5.5.d. 1460 sigma ice water mixing ratio contours and 

corresponding IR image valid at 26/12Z. Labels have been omitted from 
the ice water mixing ratio contours. 
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Figure 5.6.a. Relative humidity (with respeet to iee) vs composite ice 
water mixing ratio for 26/OOZ. Green shades indicate relative humidity 
of 90% and greater. 



Figure 5.6.b. Relative humidity (with respect to ice) vs composite ice 
water mixing ratio for 26/06Z. Green shades indicate relative humidity 
of 90% and greater. 
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Figure 5.6.C. Relative humidity (with respect to ice) vs composite ice 
water mixing ratio for 26/12Z. Green shades indicate relative humidity 
of 90% and greater. 
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VI. DISCUSSION 


Case I involved a fairly well-organized upper-level low moving into the 
southwestern United States, with an initial moisture tap provided by tropieal eonveetive 
aetivity. The ease was ehosen to study the development/evolution of the cloud ice water 
field in a region of broad upward vertical motion. Although tropical convection was 
initially responsible for injecting moisture into the upper levels, the sustainment and 
development of ice water later in the period was due to this large-scale upward motion 
ahead of the upper-level low off the coast. 

Case II focused on the model’s handling of cirrus associated with an upper-level 
ridge downstream of an upper trof. Cirrus were present in two distinct locations: 
downstream of the ridge axis (in a typically non-generative area) and downstream of the 
upper trof to the west of the ridge. The cirrus near the low was much lower in altitude 
than the cirrus downstream of the ridge axis. 

Both cases are very similar with regard to basic features, but evolution of the 
cirrus field in each case differs significantly due to ridge amplitude and available 
moisture at the source. The first case showed weaker ridge amplitude, but had a very 
strong initial moisture supply from below due to strong dynamic lifting associated with 
convection over the tropics. The second case had a very sharp ridge but the moisture 
source was not quite as strong as for the first case. 

Despite these differences, the dynamics associated with the two cases are similar 
in some respects. In the first case, tropical convective activity initially injected large 
amounts of moisture into the upper atmosphere. That moisture was carried by upper- 
level winds in advance of a closed low into an area of widespread difluence ahead of the 
low. This acted to provide additional large-scale upward vertical motion of sufficient 
strength to promote widespread nucleation before confluence downstream suppressed 
cirrus formation. In the second case, lift associated with the upper trof resulted in cirrus 
formation, and the cirrus clouds were subsequently advected across the ridge axis and 
downstream by upper-level winds before dissipating in the downstream portion of the 
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ridge. In the seeond ease, the model’s tendeney to link upward motions to iee produetion 
beeame apparent. Further study will help to determine if this is more generally the ease. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The results show that vertieal motions in the AFWA MM5 model appear strongly 
tied to presenee of iee cloud, and that the ice forecast will live or die by correct 
forecasting of vertical motions together with moisture. As we noted in the first case, 
cirrus over the Califomia/Arizona border in the satellite image shown in Figure 4.1.a 
were not forecast by the model (as shown in Figure 4.3.a). Figure 4.7.a does show 
relative humidity levels very slightly in excess of 90% over the location, so we can draw 
the conclusion that either lift was insufficient to produce ice or some microphysical 
aspect of the parameterization prevented cirrus formation. Indeed, vertical velocities 
diminished over the area through the case study period. In the second case, areas of 
cirrus were overforecast over parts of Colorado, Wyoming, Utah, and Idaho due to the 
forecast presence of upward motion and moisture in those areas. With the case studies 
presented in this thesis, it appears as though the 90 percent relative humidity field (with 
respect to ice) may be very useful as a first guess for anticipated cirrus coverage 
assuming the model is accurately handling timing and orientation of weather systems. 

This study has examined AFWA MM5 strengths and weaknesses in forecasting 
cirrus development and evolution between the 18-30 hour forecast periods for two unique 
cases of interest. As we have seen, fully explaining development and evolution of cirrus 
based on dynamic factors alone is not possible, as there are only a limited number of 
dynamic factors involved in cirrus production that can be readily analyzed. Most factors 
involved in cirrus production take place on a scale smaller than we can analyze using 
standard analysis tools such as those used in preparation of this thesis. For a better 
understanding of potential weaknesses and strengths in cirrus forecasting with the MM5, 
it is necessary to analyze many different cases and compile results for comparison. 

Perhaps the most thorough method for studying modeled production of cirrus 
would provide a researcher the ability to see the results of each term in the Reisner 
parameterization as the model ran over time, so each member’s contribution to ice water 
production could be noted. This process would be performed at each model gridpoint in 
the domain, and then the applicable dynamics would be brought into consideration along 
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with the results of the Reisner seheme. The proeess would require automation due to the 
shear volume of data involved. Onee these results were generated, the eloud mask results 
eould be overlaid with an ice water composite in a similar manner to what was done in 
this thesis, and statistics could be generated to determine errors in coverage at each 
gridpoint on the model domain. Errors could then be studied and categorized according 
to their cause, whether related to the model’s movement of features, errors in model 
dynamics, or insufficiencies in the parameterization. 

The results of this study are valid for only two cases and therefore make it 
difficult to comment on the statistical accuracy of the MM5 when predicting cirrus over a 
variety of synoptic situations. Yet, it can be said that for the two cases analyzed there 
appears to be a link between upward motions in the model, moisture in the model 
(specifically, relative humidity with respect to ice in excess of 90%), and the presence of 
model-forecast ice cloud. Additional study at a deeper level is necessary to gain an 
understanding of the effects of microphysical details on AFWA MM5 ice cloud forecasts. 
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